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1.  Motivation 


This  thesis  investigates  the  potential  for  impact  ionization  engineering  (I"E)  in  the 
nanoscale  utilizing  a  three  dimensional  (3D)  electric  field  confined  within  an  array  of  III-V 
nanopillars.  EE  has  been  exploited  in  planar  avalanche  photodetectors  with  impressive 
reductions  in  the  excess  noise,  due  to  the  size  of  the  multiplication  region  being  on  the  order  of 
the  ionization  path  length  of  the  carrier.  This  concept  can  be  extended  to  3D  nanoscale 
multiplication  regions  (3D-NMR)  where  the  injection  of  the  photogenerated  carrier  can  be 
specifically  engineered  to 


3D-APD 


impact  ionize  with  narrow 
probability  distribution 
functions  of  the  ionization 
path  lengths.  Nanoscale 
engineering  of  impact 
ionization  can  allow  an 


Nanopillar 
Optical  Antenna 
(NOA) 

X  3D  Nanoscale 
Multiplication 

Region  (so-  p|anar  APD 


unprecedented  reduction  in  Distance  from  injection  [pm] 

excess  noise  factors  and  Figure  l:  Comparison  of  engineering  the  avalanche  statistics  in  3D  NMR  compared  to 

a  planar  APD 

improvements  in  gain 

bandwidth  products.  Figure  1  illustrates  the  expected  improvement  in  avalanche  statistics  based 
on  the  proposed  3D-N0AAD).  In  a  3D-N0AAD,  photogenerated  carriers  absorbed  due  to  the 
Nanopillar  Optical  Antenna  are  precisely  swept  into  the  3D-NMR  where  they  impact  ionize 
immediately  due  to  sharp  gradients  of  the  electric  fields.  Due  to  the  geometry  of  the  3D-NMR, 


the  probability  distribution  of  the  ionization  path  length  becomes  narrower  producing  a  less 
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noisy  avalanche  gain.  This  is  in  contrast  to  a  planar  APD  where  a  planar  absorption  and  volume 
results  in  broader  ionization  path  lengths. 

1.1  Background  on  avalanche  photodetectors 

Avalanche  photodetectors  (APDs)  are  critical  components  for  a  wide  range  of 
technologies  such  as  single  photon  detectors  for  quantum  cryptography,  high  bit-rate  telecom 
receivers,  and  active  pixel  elements  for  3D  LADAR  imagers  '  .  APDs  absorb  light  generating  a 
current,  this  current  is  internally  amplified  due  a  process  called  impact  ionization.  Typically, 
APDs  are  benchmarked  in  terms  of  their  dark  current,  excess  noise  factor  (F)  and  gain- 
bandwidth  product4.  The  dark  current  is  due  to  carriers  other  than  photogenerated  carriers  being 
produced  due  to  thermal  generation,  diffusion,  and  tunneling  processes.  Excess  noise  factor  is  the 
ratio  of  the  measured  noise  spectral  density  to  the  multiplied  shot  noise5  or  can  be  understood  as 
the  fluctuations  in  the  gain.  Avalanche  build-up  time  is  the  duration  of  the  impulse  response  of 
an  APD  in  the  presence  of  gain  and  detennines  the  limit  of  the  gain-bandwidth  product.  These 
device  parameters  are  critically  dependent  on  the  electric  field  engineered  in  the  APD  and  the 
electron  and  hole  ionization  coefficients  a  and  [3  respectively.  McIntyre  predicted  that  the  excess 

noise  factor  can  be  related  to  the  mean  gain  by6  —  ,  where 

keff  =  a  /  (5  .  McIntyre's  result  predicted  that  the  excess  noise  depends  only  on  the  mean  gain  and 

keff  ,  but  not  the  size  of  the  multiplication  region.  However,  shrinking  the  volume  of  the 

avalanche  region "  and  engineering  sharp  gradients  in  the  electric  field  lowers  kefj  resulting 

in  lower  excess  noise  factor  and  high  gain-bandwidth  products4  beyond  what  is  predicted  by 
McIntyre's  model.  For  impact  ionization  to  take  place  a  carrier  must  travel  a  certain  distance 
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known  as  the  dead  space  to  gain  sufficient  energy  from  the  electric  field.  Once  the  multiplication 
volume  becomes  on  the  order  of  the  dead  space  length,  kefl  is  significantly  reduced  resulting  in  a 

lower  F  for  a  given  M. 

An  additional  benefit  of  a  lower  ke(j  would  be  a  higher  gain-bandwidth  products4. 

Emmons  predicted  the  gain-bandwidth  product  achieved  by  APDs  would  be  limited  by  the  ratio 
of  the  electron  and  hole  ionization  coefficients.  This  is  because  the  time  required  to  remove  the 
total  charge  generated  from  the  avalanche  process  increases  as  the  initial  carrier  produces  more 
secondary  carriers.  The  key  point  being  that  if  an  electron  is  initiating  the  avalanche  process, 
secondary  holes  from  the  initial  avalanche  process  could  initiate  further  electron  avalanches 
which  add  to  the  transit  time  delay.  In  traditional  planar  APDs  where  the  electric  field  is 
uniform  in  the  plane  parallel  to  the  plane  of  epitaxy,  and  the  i-region  serves  as  the  avalanche 
region  for  impact  ionization,  low  values  of  ke[j  usually  comes  at  the  cost  of  reduced  ionization 

coefficients  and  higher  breakdown  voltages.  Novel  multiplication  region  designs  are  necessary  to 
overcome  these  traditional  tradeoffs. 

The  improvements  to  APD  performance  mainly  results  from  a  more  localized  and 
coherent  impact  ionization  process  due  to  reduced  phonon  scattering.  Once  carriers  are  injected 
into  the  avalanche  region  it  is  favorable  for  them  to  ionize  with  a  specific  energy  and  distance. 
However,  phonon  scattering  broadens  the  probability  distribution  of  the  ionization  distance 
resulting  in  greater  excess  noise.  Shrinking  the  volume  of  the  avalanche  region  results  in  less 
collisions  of  these  energetic  carriers  with  phonons,  producing  a  more  coherent  impact  ionization 
process. 
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In  a  3D-N0AAD  there  is  unprecedented  flexibility  in  engineering  the  avalanche  statistics 
as  the  electric  field  distribution  is  not  only  detennined  by  the  3D  doping  profile  of  the 
semiconductor  nanostructure,  but  also  its  dielectric  environment  .  Further  investigation  of  the 
fundamental  physics  of  lower  excess  noise  in  nanoscale  multiplication  regions  is  necessary  to 
understand  the  tradeoffs  between  lower  keff  and  dark  current  due  to  tunneling. 

1.2  State  of  the  art  APD  technology 

Silicon-based  APDs  have  high  QE  in  the  visible,  however  their  efficiency  drops  significantly  in 
the  near  infrared  (IR)  wavelengths  due  to  reduced  absorption  coefficient.  Germanium  is  not 
suitable  for  IR  APDs  because  they  require  cryogenic  cooling  (77K)  to  reduce  the  high  dark 
current.  Thus,  InGaAs/InP  APDs  are  the  workhorse  for  near  infrared  photodetection  with  QE's  of 
~  90%  and  excess  noise  factors  of  ~  3  at  a  gain  M  =10  '  . 

Regarding  application,  the  operating  reverse  bias  in  state-of-the-art  InP-InGaAs  APDs  is 
typically  in  excess  of  40  V.  This  value  is  just  slightly  less  than  the  large  breakdown  voltage 
required  to  maintain  acceptable  levels  of  field-assisted  dark  carriers  in  the  long  multiplication 
region  (>1.5  pm).  Lowering  the  reverse  bias  to  below  -  15  V  is  desirable  as  it  simplifies  the  bias 
circuits,  especially  in  array  configurations.  To  reduce  the  breakdown  voltage  in  a  planar  APD 
requires  the  multiplication  region  to  be  below  500  mn.  However,  the  required  electric  field  to 
meet  the  breakdown  condition  increases  field-assisted  dark-carrier  generation  to  unacceptable 
levels.  Thus,  wide  multiplication  layers  (>1.5  pm)  are  preferred  in  planar  APDs  to  minimize 
dark  counts  at  the  expense  of  the  large  reverse  biases. 
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1.3  Active  imaging  and  photon  counting  LADAR  technology 

One  of  the  main  application  areas  of  3D  NOAADs  are  LADAR  imagers,  which  allow  the  photon 
arrival  times  from  objects  to  be  determined  allowing  a  3D  correlation  of  the  object  space  10  - 
100  km  away5, 12.  Fig.  2a)  below  shows  a  schematic  of  3D  imaging  system16  where  a  reflective 


30 


Figure  2:  Schematic  of  3D  imaging  system  (taken  from  Albota  et.  al16)  b)  Probability  of  single  photon  detection  vs.  receiver  noise 
threshold  plotted  for  different  values  of  kef/2 

cone  is  illuminated  by  pulsed  laser  illuminator.  The  reflected  photons  are  detected  are  sensed  by 
an  array  of  single  photon  detectors  where  the  time  of  flight  instead  of  the  intensity  of  photons 
allows  the  creation  of  a  3D  image  of  the  cone.  The  number  of  photons  arriving  at  the  APD  array 
is  extremely  small  such  that  the  APD  array  requires  sensitivity  to  single  photon  counting  levels. 
In  the  case  of  a  linear  mode  single  photon  detector,  keff  has  a  very  large  impact  on  the  probability 
of  detecting  a  single  photon.  From  Fig.  2b)  it  can  be  seen  that  the  probability  of  detecting  a 
single  photon  is  ~  60%  given  a  keff  ~  0  and  discrimator  threshold  of  60.  Single  carrier 
multiplication  devices  are  therefore  extremely  useful  for  single  photon  detection  with  high 
detection  probability. 
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A  LADAR  imager  consists  of  APD  arrays  which  are  either  operated  in  Geiger  mode  or  linear 
mode.  Geiger  mode  involves  biasing  the  APD  above  its  breakdown  voltage  such  that  the  output 
current  saturates  while  linear  mode  involves  biasing  the  APD  below  the  breakdown  voltage  such 
that  output  current  is  linearly  proportional  to  the  primary  photocurrent.  Geiger  mode  LADAR 
imagers  have  been  preferred  for  single  photon  sensitivity,  however  they  come  at  the  cost  of 
spurious  detection  events  due  to  dark  carriers.  These  spurious  detection  events  slows  down 
Geiger  mode  imagers  limiting  detailed  3D  imaging.  Linear  mode  operation  on  the  other  hand 
don’t  suffer  from  spurious  detection  events  but  needs  to  be  implemented  in  concert  with  high 
perfonnance  amplifiers  to  reach  single  photon  detection  sensitivity. 

LADAR  imagers  work  on  the  principle  of  a  pulse  laser  illuminator  generating  laser 
pulses  which  reflect  off  a  target  and  are  subsequently  detected  by  the  APDs.  The  returned  pulses 
are  imaged  onto  each  pixel  to  form  an  image  and  sampled  to  obtained  3D  infonnation  of  the 
object.  LADAR  imagers  based  on  Geiger  mode  APDs  have  been  demonstrated  .  The  main 
problem  with  state  of  the  art  Geiger  mode  APDs  is  that  single  photon  detection  comes  at  the 
price  of  large  dark  currents  and  long  deadtimes.  Sensitivity,  speed  and  bandwidth  requirements 
of  the  APD  drive  the  design  of  LADAR  imagers. 

1.4  Subwavelength  metallic  antennas  coupled  nanophotodiodes 

Reduced  quantum  efficiency  is  a  major  challenge  faced  when  scaling  down  the 

photodetector  volume.  A  possible  solution  comes  from  the  radio -frequency  (RF)  frequency 

range,  where  metallic  antennas  in  receivers  concentrate  the  electromagnetic  field  within  a 

subwavelength  volume.  Researchers  have  explored  different  antenna  configurations  in  coupling 

light  into  subwavelength  absorption  volumes.  Rectangular  gold  nanorod  antenna  array  on  silicon 
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18,  circular  metal  hole  arrays  coupled  to  intersubband  transitions  in  InAs  QDs19  and  free  standing 
metallic  gratings  integrated  with  HgCdTe  photodetectors  exemplify  the  utilization  of  extra 
ordinary  optical  transmission  (EOT)  using  near  field  optical  field  enhancements. 

Indications  that  these  concepts  could  be  applied  to  the  infrared  and  optical  regime  were  further 
encouraged  by  observations  of  so-called  “extraordinary  optical  transmission”  of  light  through 
arrays  of  sub-wavelength  holes  in  a  metal  film  21'27,  or  a  single  sub-wavelength  aperture 
surrounded  by  a  periodic  bullseye  pattern  “  '  .  This  phenomenon  is  accompanied  by  a  significant 
enhancement  of  the  electric  field  within  the  subwavelength  aperture,  which  can  be  leveraged  for 

25  32  33 

a  variety  of  photonic  applications  ’  ’  . 

In  general,  the  benefit  of  plasmonically-enhanced  photodetection  is  that  the  “plasmonic  antenna” 

’  '  ,  effectively  collects  light  from  a  large  area  and  concentrates  it  on  a  much  smaller  and 

thinner  photodiode  junction  resulting  in  lower  capacitance,  and  lower  transit  time.  Various  field 

enhancement  mechanisms  can  be  leveraged,  including  the  enhancement  near  the  metal/dielectric 

interface  in  a  propagating  Surface  Plasmon  Polariton  Bloch  wave  (SPP-BW)  36’  40.  Also, 

resonant-field  enhancement  occurs  in  bounded  metallic/dielectric  structures  that  support 

Localized  Surface  Plasmon  Resonances  (LSPRs)34’  35,  such  as  metallic  nanoparticles  or 

nanoholes.  These  resonant  field  “focusing”  mechanisms  effectively  increase  the  optical  path 

length  to  be  many  times  larger  than  the  physical  junction  thickness.  In  principle,  this  allows  the 

reduction  of  the  detector  volume  even  below  the  diffraction  limit  for  more  efficient  integration, 

reduced  leakage  current  for  greater  detector  sensitivity,  and  reduction  of  the  junction  capacitance 

and  transit  time  for  higher  speed  detectors.  Again,  these  demonstrations  are  planar  (2D) 

configurations  and  based  on  complex  lithographic  techniques.  We  note  that  planar  demonstration 
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of  the  plasmonic  detectors  fall  substantially  below  theoretical  predictions  due  to  overlap  of  field 
intensity  “hot-spots”  with  lossy  metal.  However,  detectors  operating  in  the  near  infrared  ( — 1 .55 
pm)  are  spectrally  close  to  the  plasma  resonances  of  noble  metals  and  are  expected  to  have  large 
metal  losses.  This  is  not  a  concern  in  the  case  of  the  NOA  due  to  the  tight  focusing  of  light 
within  the  semiconductor  nanopillar  instead  of  the  lossy  metal.41 


1.4  Promise  of  proposed  3D-NOAADs  -  relaxing  the  QE  vs.  dark  current  tradeoff 

In  comparison  to  planar  APDs,  the  structure  and  electric-field  distribution  of  the  3D-APD 
have  several  key  advantages  with  regards  to  impact  ionization  and  APD  sensitivity  enhancement. 
The  most  important  figure  of  merit  to  be  improved  in  an  APD  is  the  detectivity  which  quantifies 


the  improvement  in  sensitivity  as  defined  in  Equation  1  as 


42 


D  =■ 


A  77 


2hc(G(T)t ) 


f  ,  A 


(1) 


where  h  is  the  Plank  constant,  c  is  the  speed  of  light,  A  is  the  wavelength  of  light,  77  in  the 
quantum  efficiency,  Aa  is  the  optical  collection  area,  Ae  is  the  electrical  area  of  the 
photodetector  absorber,  t  is  the  thickness  of  the  photodetector  absorber,  G(T)  is  the  thermal 
generation  rate.  From  equation  1  it  can  be  seen  that  the  specific  detectivity  of  the  photodetector 
can  be  improved  by  1)  decoupling  the  optical  collection  from  the  electrical  area  such  that  Ao 

»  Ae  and  2)  decreasing  the  volume  of  the  photodetector  absorption  material  such  that  the  due 


generation  noise  ( Gt )2  can  be  reduced.  The  3D-NOAD  focuses  light  from  a  large  collection 
area  (  Ao )  to  a  small  absorbing  foot-print  (  Ae )  within  the  nanopillar  enabling  high  QE.  The  small 
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footprint  enables  nanoscale  engineering  of  the  impact  ionization  process  for  improved  APD 
figures  of  merit. 
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2.0  Nanopillar  Optical  Antenna  for  efficient  wavelength  tunable 
absorption 

The  NOA  is  the  crucial  technology  enabling  footprint  reduction  of  the  3D-NMR,  while  maintaining  high 
quantum  efficiency.  The  NOA  is  a  3D  plasmonic  antenna  structure  which  collects  electromagnetic  energy 
and  concentrates  it  in  a  subwavelength  area  by  coupling  with  surface  plasmon  waves21'27' 4j.  In  the  optical 
frequency  range  where  the  dielectric  function  of  metals  take  on  the  plasma  form  these  surface  waves  are 
known  as  Surface  Plasmon  Polaritons  Bloch  Wave  (SPP-BW)21’  44  or  Localized  Surface  Plasmon 
Resonances  (LSPRs).  SPP-BWs  have  in-plane  momentum,  and  therefore  needs  to  be  phase-matched  by  a 


Figure  3:  a)  The  Nanopillar  Optical  Antenna  allows  the  coupling  of  both  SPPs  and  LSPRs  via  geometry  partial  metal  shell  and  self 
aligned  metal  nanohole  b)  Calculated  external  quantum  efficiency  of  the  NOA  as  a  function  of  pitch  of  the  self -aligned  nanohole 
lattice 


periodic  surface  corrugation  to  couple  to  normally  incident  light.  LSPRs  are  supported  by  the  specific 
shape  of  the  metallic  nanostructure22’ 45  and  can  directly  couple  to  normally  incident  light.  This  thesis  will 
demonstrate  that  the  NOA  can  utilize  both  SPP-BWs  and  LSPRs  for  harvesting  light  for  enhanced 
absorption  within  the  nanopillar41’ 46. 
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The  NO  A  consists  of  a  partial  gold  shell  on  the  nanopillar  sidewalls  and  gold  cap  on  top  of  the  nanopillar 
along  with  a  self  aligned  nanohole  as  shown  in  Fig.  2  a).  The  width  of  the  nanopillar  (w)  and  height  (h) 
allows  the  focussing  of  LSPRs  while  the  periodicity  (P)  of  self  aligned  nanohole  allows  focussing  of  SPP- 
BWs  within  the  nanopillar  as  shown  in  Fig.  2  b).  The  optical  field  enhancement  due  to  the  NOA  takes 
place  in  the  nanopillar  instead  of  the  lossy  metal  (as  shown  in  Fig.  2  b)),  which  is  the  key  enabling 
technology  for  nanoscale  footprint  and  high  quantum  efficiency  of  the  3D-NOAAD.  The  geometry  of  the 


Figure  4:  a)  FDTD  simulations  of  the  power  absorption  spectra  b)  Measured  photoresponse  spectra 

NOA  can  be  tuned  for  high  external  quantum  efficiencies  at  1  pm.  Figure  1  c)  shows  a  contour  map  of  the 
quantum  efficiency  of  the  NOA  as  a  function  of  periodicity  (P)  while  w  and  h  are  fixed.  The  NPs  consist 
of  lno.35Gao.47As  with  a  diameter  of  225  nm,  and  550  nm  in  height  arranged  in  square  lattice.  The  dotted 
white  and  black  lines  corresponds  to  first  order  and  second  order  SPP-BWs.  In  addition  there  is  a  third 
mode  around  X  =  1  pm  which  is  independent  of  pitch,  this  is  a  LSPR  excited  due  to  the  width  and  height 
of  the  nanopillar. 

Due  to  the  hybridization  of  the  second  order  SPP-BW  mode  and  LSPR  a  peak  quantum  of  ~  50%  is 
achieved.  These  initial  calculations  indicate  that  the  NOA  can  produce  an  external  quantum  efficiency 
in  excess  of  ~50%  shown  in  Fig.  3a)  at  1  pm.  Increasing  the  pitch  of  the  nanopillar  results  in  the  peak 
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absoiption  of  the  resonances  to  shift  as  shown  in  Fig.  3b).  The  ability  to  integrate  tunable  light  absorption 
using  the  NOA  opens  applications  for  hyperspectral  imaging  in  SW1R  and  MW1R  wavelengths. 

Reflection  spectroscopy  of  plasmonic  antenna  resonances 

A  key  aspect  of  the  electromagnetic  design  of  the  NOA  structure  is  to  critically  couple  incident  radiation 
to  the  semiconductor  nanopillars.  In  order  to  characterize  optical  absoiption,  simultaneous 
reflection/transmission  measurements  need  to  be  performed.  Typically  reflection/transmission 
measurements  are  performed  using  unfocussed  collimated  light  on  large  areas  on  the  order  of  ~  cm2. 
However,  since  NOA  is  manufactured  using  electron  beam  lithography  it  is  not  feasible  to  produce  NOA 
arrays  in  a  large  cm2  area. 

The  most  significant  challenge  in  measuring  optical  absoiption  is  focusing  the  excitation  light  to  the  NOA 
patterned  area  on  the  order  of  50  pm  x  50  pm.  Tightly  focusing  incident  light  using  high  numerical 
aperture  (NA)  objectives  generates  light  with  a  range  of  incident  angles.  For  example,  typical  Swazichild 
reflective  objectives  with  a  NA  of  0.65  have  incident  angles  ranging  from  17°  -  42°  broadening  plasmonic 
absoiption  features.  Despite  the  broadening  of  plasmonic  resonances,  focused  reflection  measurements 
offer  the  highest  spatial  resolution  characterization  limited  by  the  spot  size47  of  the  source  or  aperture  size 
of  the  microscope48. 

In  the  Short  Wave  Infrared  (SW1R)  there  are  two  schemes  for  focused  reflection  characterization:  1)  A 

Fourier  transform  infrared  (FT1R)  spectrometer  (Thermo  Scientific  8700)  integrated  with  an  1R 

microscope  (Thermo  Scientific  continuum)  and  a  2)  fiber  coupled  spectrometer  with  NIR  objective  such 

as  filmetrics  F20-EXR.  In  both  cases  a  white  light  source  allows  spectral  measurements  in  the  near 

infrared  (NIR).  The  main  distinction  between  the  two  systems  is  that  the  fiber  coupled  filmetrics  system  is 

limited  to  reflection  measurements,  while  the  Thermo  Scientific  system  can  measure  both  reflection  and 

transmission  from  1-10  pm  wavelengths.  While  measurements  of  absoiption  will  give  information  only 

about  coupling  to  SPP-BW  and  LSPR  modes  regardless  of  their  coupling  to  the  semiconductor  nanopillar, 
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simultaneous  measurement  and  comparison  of  the  photocurrent  spectrum  will  provide  a  method  to 
identify  which  SPP-BW  modes  of  the  plasmonic  crystal  result  in  the  maximum  photogeneration  within 
the  nanopillar.  Fig.  4a)  and  b)  shows  a  schematic  of  the  FTIR/IR  microscope  configured  in  the  reflection 


Figure  5:  a)  Nicolet  8700  FTIR  pectrometer  combined  with  continuum  IR  microscope  b)  Reflection  mode  configuration  of  continuum  IR 
microscope  (taken  from  Thermo  Scientific  brochure)48  b)  Reflection  measurement  of  NOA  array 

mode  for  focused  measurements49.  Light  from  the  white  light  source  is  focused  onto  the  sample  using  one 
half  of  the  aperture  and  the  reflected  light  is  collected  and  sent  back  to  the  spectrometer  for  spectral 
analysis.  The  spot  can  be  observed  via  the  beam  splitter  on  the  sample.  Reflectance  measurements  were 
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performed  on  the  NOA  arrays  as  shown  in  Fig.  4c).  Both  SPP-BW  and  LSPRs  can  be  observed  as 
reflection  dips.  The  LSPR  remains  at  fixed  spectral  position  of  ~  1000  nm  as  the  pitch  is  increased  while 
the  reflection  dip  due  to  the  SPP-BW  mode  shifts  by  ~  50  nm. 


Electromagnetic  simulations  of  surface  plasmon  resonances  in  nanopillars 

To  design  the  optical  absorption  in  Nanopillar  Optical  Antenna  Avalanche  Detectors  full 
wave  finite  difference  time  domain  (FDTD)  simulations  are  carried  out  using  the  commercial 
software  Lumerical.  The  incident  plane  wave  is  polarized  in  both  the  x-  and  y-  polarizations  over 


Wavelength  [nm]  Wavelength  [nm] 


Figure  6:  Structure  definition  for  3D  FDTD  simulation,  b)  FDTD  simulation  of  the  b)  power  absorbed  in  the  x-  and  y-  polarizations,  and  c) 
Electric  field  intensity. 

a  wavelength  range  of  900  -  1600  nm.  The  3D  unit  cell  of  the  modeled  structure  is  shown  in 
figure  4  a).  Periodic  boundary  conditions  are  implemented  in  the  x-  and  y-  boundaries  with  a 
pitch  of  1  pm.  Perfectly  Matched  Layer  (PML)  boundary  conditions  are  implemented  in  the  z- 
directions.  Ino.4Gao.6As  optical  parameters50  are  used  to  describe  the  NP,  GaAs  optical 
parameters51  are  used  for  the  substrate.  A  fine  mesh  of  1  nm  x  1  nm  x  2  nm  is  used  in  the  gold 
and  the  tip  of  the  exposed  nanopillar  in  contact  with  the  gold. 

The  fractional  power  absorbed  in  the  nanopillar  is  computed  by  numerically  evaluating 
Equation  (2),  where  s"(co)  is  the  imaginary  part  of  the  dielectric  constant  of  Ino.4Gao.6As,  and 
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\e(g>,x,  y,z)|  is  the  electric  field  intensity.  The  integral  is  evaluated  within  the  volume  of  the 
Ino.4Gao.6As  nanopillar. 

pabs  (®) = \\\\co\E{co’x’y’z)\  s"(coyv  (2) 

Figure  4b)  shows  the  spectra  for  the  fractional  power  absorbed  in  the  NP  for  the  x-  and  y- 
polarizations  respectively.  For  the  x-  polarization  there  is  a  peak  at  1030  nm  where  47%  of  the 
incident  light  is  absorbed  in  the  nanopillar,  while  for  the  y-  polarization  there  are  two  peaks  at 
1010  nm  and  1060  nm.  The  peaks  at  1030  nm  for  the  x-polarization  and  1010  for  the  y- 
polarization  are  attributed  to  SP{±]  ±1|  mode.  The  SPf±]  +lj  mode  is  excited  in  both  the  x-  and  y- 

polarizations,  where  it  is  stronger  in  the  x-  than  the  y-  polarization. 

Figure  4c)  shows  the  spectra  for  the  electric  field  intensity  integrated  over  the  NP  volume 
for  the  x-  and  y-  polarizations  respectively.  Mirroring  the  calculated  power  absorption  spectrum 
in  Fig.  4b),  the  SFj±1  ±l)  are  be  observed  in  the  integrated  electric  field  intensity  for  both  the  x- 
and  y-  polarizations.  In  addition  the  x-  polarization  has  an  even  stronger  resonance  at  1500  nm 
which  is  attributed  to  the  resonance,  which  is  largely  absent  for  the  y-polarization. 

However,  this  peak  is  not  observed  in  the  absorption  spectra,  as  it  lies  at  a  longer  wavelength 
than  the  cut-off  of  the  InGaAs  NPs  at  1450nm.  The  strong  SP{U  { (  resonance  is  very  encouraging 

for  future  1.55  pm  device  demonstrations  since  greater  optical  coupling  efficiency  of  the  SP[()±y) 
mode  into  the  NP  would  result  in  a  much  higher  responsivity. 
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3.  Fabrication  of  Nanopillar  Optical  Antenna  Avalanche  Detectors 

For  the  devices  studied  in  this  dissertation,  Zn-doped  Ino.4Gao.6As  NPs  are  grown  by 
selective-area  epitaxy  (SAE)  on  n+  doped  GaAs  (111)  B  substrates  at  10  cm'  .  Figure  5  (a) 
shows  a  typical  NP  array  before  processing,  and  Fig.  5(b)  shows  an  SEM  image  of  the  fully 
processed  photodiode  array.  A  schematic  of  the  final  device  is  shown  in  Fig.  5(c).  The  NPs  have 


Figure  7:  a)  Ino.4Gao.6As  NPs  grown  on  n+  GaAs  by  selective  area  MOCVD  arranged  in  a  square  lattice  with  a  l(im  pitch,  b)  Final 
fabricated  APD  showing  high  speed  contacts  c)  Close  up  view  of  the  nanopillar  array  showing  the  signal  contact  aligned  to  the 
nanopillar  array. 

a  height  of  1.3  pm,  a  diameter  of  200  nm,  and  are  arranged  in  a  1  pm  pitch  square  lattice.  The 
NP  doping  concentration  is  estimated  to  be  5x10  cm"  from  single  pillar  resistivity 
measurements.  A  more  detailed  description  of  substrate  patterning  and  NP  epitaxy  are  published 
elsewhere52. 

Following  epitaxy,  the  NP  array  is  planarized  using  Benzocyclobutene  (BCB),  which  is 
hard-cured  and  etched  back  to  expose  -500  nm  of  the  NP  tips.  The  top  metal  contact  (Cr:  20  nm, 
Au:  200  nm)  is  evaporated  at  35  0  from  the  substrate  normal  to  coat  both  top  and  the  exposed 
side  of  each  NP  along  with  the  entire  surface  interconnecting  the  pillars.  The  NP  “shadow” 
leaves  the  other  pillar  side  uncoated  and  forms  a  self-aligned  nanohole  adjacent  to  each  NP.  The 


nanohole  length  depends  on  the  exposed  tip  length  and  the  incident  angle  of  metal  deposition, 
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while  the  nanohole  width  is  defined  by  the  NP  diameter.  The  resulting  nanoholes  are  260  nm 
wide  and  315  nm  long.  This  metal  nanohole  array  forms  a  3D  Nanopillar  Optical  Antenna  self 
aligned  to  the  NPs  described  in  Section  2. 

3.1  Polarization  effects  of  plasmonically  enhanced  nanopillar  photodetectors 

Effects  of  the  asymmetric  self-aligned  nanohole  are  elucidated  through  polarization- 
specific  spectral  response  of  the  NOAAD  described  in  Fig.  6  (a)-(c).  As  shown  in  Fig.  6(a)  the 
polarization  of  light  with  respect  to  the  long  axis  of  the  nanohole  is  considered  x-polarized  and 
the  polarization  orientation  is  delineated  using  angle  (f>. 
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Figure  8:a)  High  resolution  SEM  of  the  nanopillar  in  the  nanohole  and  orientation  of  the  light  polarization  b)  Photoresponse  for  x- 
and  y-  polarized  light  c)  Polarization  dependent  photocurrent  at  1100  nm  showing  cos2(ip) 


As  shown  in  Fig.  6(b)  the  photoresponse  of  the  device  when  incident  light  is  x-polarized 
light  (<j)  =  0  °)  produces  a  higher  photocurrent  compared  to  the  y-polarized  light  (tj)  =  90  °)  for 
much  of  the  photoresponse  spectrum,  except  at  1260  nm  where  the  photocurrent  for  both 
polarizations  is  the  same.  It  is  interesting  to  note  that  the  peak  responsivity  for  both  (|>  =  0  0  and  (j> 
=  90  0  occur  at  1110  nm  independent  of  polarization  suggesting  photoresponse  enhancement  due 
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to  SP^±l  =|)  mode.  Fig.  6c)  shows  a  polar  plot  of  the  photocurrent  for  incident  wavelength  at  1100 

nm.  The  photocurrent  is  measured  at  the  peak  responsivity  of  X  =  1100  nm  from  cj)  =  0  °  through 
a  full  360  0  revolution  at  10  0  increments.  The  maximum  photocurrent  is  measured  for  the  x- 
polarization  with  a  25%  drop  in  the  photocurrent  for  y-polarized  light.  The  photocurrent  shows  a 
cos2  (p  dependence  on  the  incident  light  polarization  similar  to  the  optical  antenna  behavior 
observed  with  gold  nanorods  .  In  the  case  of  the  control  device  there  is  no  dependence  of  the 
photoresponse  on  polarization.  This  dipole-like  dependence  of  the  photocurrent  is  commonly 
attributed  to  the  specific  nature  of  the  Localized  Surface  Plasmon  (LSP)  that  is  excited  in  the 
gold  nanohole22,  44'  53 .  However,  the  polarization  dependence  of  the  spectral  response  is  the 
opposite  expected  of  an  elongated  nanohole,  where  the  maximum  transmission  occurs  at  <j)  =  90 
and  the  minimum  at  <j>  =  0  °. 

Localized  Surface  Plasmon  Absorption 

The  NOA  is  advantageous  for  improving  the  signal  to  noise  ratio  of  nanopillar 
photodetectors  because  it  allows  optical  absorption  greater  than  geometrical  cross-section  of  the 
nanopillar.  Fig.  7(a)  shows  a  schematic  of  a  single  NOA  which  consists  of  a  nanopillar  with 
width  w,  exposed  nanopillar  height  of  h,  and  the  partial  gold  shell.  Specifically,  h  is  the  height  of 
the  nanopillar  above  the  polymer  (bisbenzocyclobutene  -  BCB)  planarization  layer.  LSPRs  can 
be  spectrally  tuned  by  varying  the  geometry  of  the  (NOA)  in  tenns  of  w  and  h.  The  NOA 
supports  LSPRs  due  to  the  metal  cap  on  top  of  the  nanopillar,  and  the  metal  shell  on  the 
sidewalls  of  the  nanopillar.  The  absorption  of  a  single  NOA  is  studied  by  full  wave  FDTD 
simulations  using  the  total  field  scattered  field  approach.  In  the  TFSF  approach  the  boundary 
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conditions  are  Perfectly  Matched  Layer  (PML)  boundary  conditions  in  the  x-  y-  and  z- 
boundaries,  and  the  illumination  is  replaced  with  a  total  field/scattered  field  source,  which 
isolates  the  scattered  light  outside  the  source  54, 55. 

The  optical  properties  of  gold  in  the  near  infrared  was  modeled  using  Johnson  and 
Christy  56 .  Light  is  normally  incident  and  polarized  along  the  long  edge  (x-direction)  of  the 
nanohole,  such  that  it  is  more  efficiently  absorbed  as  reported  previously  56.  The  total  electric 
field  intensity  enhancement  inside  the  nanopillar  is  quantified  using  the  normalized  absorption 
cross-section  ( Qabs ),  which  is  defined  as  the  absorption  cross-section  normalized  to  the  cross- 
section  of  the  nanopillar  as  shown  in  Eq.  (4).  The  parameters  of  Eq.  (4)  are  co  which  is  the 
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Figure  9:  a)  Schematic  of  Nanopillar  Optical  Antenna  (NOA).  Contour  plot  of  the  Qabs  as  b)  w  is  varied,  and  c)  as  h  is  varied. 


angular  optical  frequency,  s"  is  the  imaginary  part  of  the  dielectric  permittivity  of  Ino.3Gao.7As50, 

I  _i2 

\E\  is  the  electric  field  intensity  inside  the  nanopillar,  Isource  the  intensity  of  the  total  field 

scattered  field  source,  and  ANP  the  cross-section  area  of  the  hexagonal  nanopillar  54.  Qabs  is  a 
measure  of  the  funneling  of  light  from  a  larger  collection  area  and  absorption  into  the  nanopillar. 
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Fig.  7(b)  and  (c)  show  contour  plots  of  Qabs  as  the  geometry  of  the  NOA  is  varied.  In  Fig. 
1(b)  h  is  kept  constant  at  180nm,  and  the  width  of  the  nanopillar  is  varied.  The  peak  value  of  Qabs 
smoothly  shifts  to  longer  wavelengths  with  increasing  width  such  that  a  maximum  value  of  -3.9  is 


900  1000  1100  1200  1300  1400 

Wavelength  [nm] 


300 


18 
15 
12 
-  9 

g  6 

3 
0 


Figure  10:  a)  Absorption  enhancement  due  to  the  Localized  Surface  Plasmon  (LSPR)  supported  by  the  metal  cap  b)  Electric  field 
intensity  profile  of  LSPR  mode  at  1270  nm  in  x-z  plane 

obtained  for  a  width  of  -  160  nm.  In  Fig.  7(c)  w  is  kept  constant  at  160  nm,  and  the  exposed  height  h 
is  varied.  The  peak  value  of  Qabs  in  this  case  also  shifts  to  longer  wavelength  as  the  exposed  height 
increases,  however  there  are  two  hot  spots  in  the  absorption  cross-section.  The  first  hot  spot  in  the  h 
=  200  -  250  nm  range  is  due  to  an  LSPR  supported  mainly  by  the  metal  cap.  The  second  hot  spot  in 
the  h  =  300  -350  nm  is  due  to  an  LSPR  supported  mainly  the  partial  metal  shell.  These  simulations 
show  that  the  NOA  supports  LSPRs,  and  that  the  optical  absorption  can  be  tuned  by  the  three 
dimensional  geometry  of  both  the  nanopillar  gold  cap  and  partial  gold  shell  covering  the  nanopillar. 


20 


34 


The  highest  enhancement  in  Qabs  is  obtained  when  h  is  equal  to  the  thickness  of  gold  on  polymer 
planarization  such  that  w  =  160  nm  and  h  =  180nm. 

A  representative  spectrum  of  the  enhancement  of  Qabs  is  shown  in  Fig.  8(a).  The  NOA 
geometry  is  such  that  w  =  160  nm  and  h  =  320  nm.  The  peak  value  of  Qa),s  is  ~  2.8,  while  the  peak 
absorption  takes  place  at  1270  nm.  This  geometry  of  the  NOA  was  chosen  such  that  the  LSPR  could 
be  clearly  distinguished  from  the  SPP-BW  modes.  Fig.  8(b)  shows  the  electric  field  intensity  in  the 
x-z  plane.  The  nanopillar  is  outline  by  a  dotted  white  line.  Electric  field  enhancements  within  the 
nanopillar  can  result  in  useful  photogeneration  and  enhanced  photocurrent.  The  electric  field  extends 
from  the  metal  cap  on  top  of  the  nanopillar  down  to  the  nanopillar  where  it  is  absorbed.  The  LSPR 
due  to  the  metal  cap  can  be  tuned  by  varying  the  exposed  height  of  the  nanopillar.  At  its  peak  value, 
the  electric  field  intensity  inside  the  nanopillar  is  xlO  the  incident  electric  field  intensity.  The  LSPR 
forms  a  vertical  dipole  resonance  from  the  metal  cap  on  top  of  the  nanopillar  and  the  partial  metal 
shell  on  the  nanopillar  sidewalls.  This  kind  of  LSPR  is  qualitatively  similar  to  nanopatch  antenna 
structure  used  for  plasmonic  lasers  except  the  case  of  the  NOA  has  a  broken  ground  plane  57. 
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Figure  11:  .  a)  Electric  field  intensity  in  the  y-z  plane  of  the  LSPR  mode  at  1270  nm  and  the  corresponding  electric  field  in  b)  the  z- 
direction  Ez  and  c)  y-direction  Ey 
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The  nature  of  the  LSPR  due  to  the  partial  gold  shell  can  be  further  investigated  by 
decomposing  the  electric  field  intensity  into  its  constituent  field  components.  Fig.  9(a)  shows  the 
electric  field  intensity  of  the  NOA  structure  in  the  y-z  plane.  The  electric  field  intensity  enhancement 
is  greatest  near  the  partial  gold  shell  on  the  nanopillar  sidewalls  (greater  than  xl8  the  incident 
electric  field  intensity),  and  lowers  towards  the  center  of  the  nanopillar.  Fig.  9(b)  shows  the  electric 
field  component  in  the  z-direction  (Ez).  Ez  is  pointing  downwards  in  the  nanopillar  and  maximum  in 
the  center  of  the  nanopillar  .  Fig.  9(c)  shows  the  electric  field  component  in  the  y-direction.  The 
direction  of  Ey  changes  on  the  side  walls  of  the  nanopillar  due  to  the  partial  metal  shell  coating  the 
sidewalls.  These  results  indicate  that  the  enhancement  in  electric  field  intensity  of  the  NOA  is  due  to 
LSPRs  of  the  metal  cap  on  the  nanopillar  producing  enhancement  in  the  z-directed  electric  field  and 
the  partial  metal  shell  on  the  sidewalls  of  the  nanopillar  producing  an  enhancement  in  the  y-directed 
electric  field. 

Angular  spectral  response  mapping  of  plasmonic  absorption 

Fig.  1 1  shows  the  experimental  angle-dependent  photoresponse  at  zero  bias.  For  the 

photoresponse  characterization  the  photocurrent  of  the  wire  bonded  device  is  measured  using  a 

Keithley  2402  multimeter.  Light  from  a  quartz  halogen  lamp  is  dispersed  by  a  grating 

monochromator  and  is  incident  on  the  device.  In  order  to  characterize  the  plasmonically 

enhanced  absorption,  photocurrent  contributions  from  the  different  plasmonic  modes  need  to  be 

identified  and  quantified.  SPP  modes  have  a  characteristic  dispersion  related  to  its  propagation  in 

the  plasmonic  crystal.  This  dispersion  can  be  measured  using  variable  angle  photocurrent 

spectroscopy,  as  illustrated  in  Figure  1 1(a)  19’ 41’ 58’ 59.  This  measurement  can  be  used  to  quantify 

what  proportion  of  incident  radiation  is  coupling  to  SPP  modes  and  map  these  modes  by  varying 
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Figure  12:  a)  Schematic  of  the  sample  orientation  in  the  angle  dependent  photoresponse  measurement,  b)  Angular  photoresponse  of 

the  gold  nanohole  array  photodetector  showing  the  shift  and  splitting  of  the  SP  resonance. 

u 

the  angle  of  incidence.  A  quartz-tungsten  halogen  lamp  is  used  as  the  illumination  source  in  a 
variable  angle  photocurrent  spectroscopy  setup  that  allows  measurements  spanning  incident 
angles  from  0-80°.  Typical  angular  photoresponse  spectra  of  plasmonically  enhanced  nanopillar 
detector  is  shown  in  Figure  10(b).  Specifically,  Figure  10(b)  shows  that  the  peak  responsivity 
occurs  at  — 1110  nm  is  roughly  described  by  the  (1,1)  scattering  order  of  the  SPP  mode  . 

SP{  i  ,  j  degenerate  modes  will  split  into  the  SP±i,.i  and  SP±ij+i  modes  with  increasing  incident 

angle  9.  At  normal  incidence  (9  =  0)  the  SP±i,  ±i  modes  are  degenerate,  resulting  in  the  peak  at 
1112  nm.  As  9  is  increased  the  photoresponse  peak  at  1112  nm  red  shifts  due  to  the  SP±i,-i 
modes,  and  also  broadens  due  to  the  SP±ij+i  modes.  Eventually  at  9  >  15  °,  the  single  peak  has 
clearly  split  into  two  peaks  in  good  qualitative  agreement  with  the  interpretation  of  the 
photoresponse  shift  as  being  due  to  SPP  field  enhancement  in  the  nanopillars.41 
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2.1  Hybridizing  LSPRs  and  SPP-BWs 

NOAs  can  be  arranged  in  a  lattice  such  that  the  optical  absorption  can  be  further 
enhanced  by  the  overlap  of  SPP-BWs  56  and  LSPRs  resulting  in  hybridized  plasmonic  modes  for 
enhanced  photoresponse.  The  boundary  conditions  of  the  simulation  are  changed  to  periodic 
boundaries  in  the  x-  y-  directions,  while  the  incident  illumination  is  a  plane  wave  source.  The 
nanopillar  exposed  height  and  nanopillar  width  are  chosen  to  be  h  =  320  mn  and  w  =  160  nm  to 
couple  LSPRs  due  to  the  NOA  structure  described  earlier.  A  pitch  of  580  nm  was  chosen  such 
that  1)  the  resonance  of  the  SPP-BW  can  be  clearly  distinguished  from  the  NOA  resonance  and 
2)  resonances  are  within  the  absorption  spectrum  of  the  Ino.35Gao.65As  material. 

Fig.  11(a)  shows  the  fractional  power  absorbed  in  the  nanopillar  per  unit  cell.  The  first 
peak  at  1000  nm  is  due  to  the  SPP-BW  mode,  which  is  only  appears  when  periodic  boundary 
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Figure  13:  a)  The  full-wave  simulation  of  the  optical  absorption  in  the  nanopillar  array  for  P  =  580  nm.  Power  absorbed  in  the  x-z 
plane  for  b)  the  SPP  mode  and  c)  LSPR  mode. 

conditions  are  used.  While  the  second  peak  at  1300  nm  is  due  to  the  LSPR,  which  is  the  same  as 
is  present  in  Fig.  8(a).  Approximately  18%  of  the  incident  light  is  absorbed  in  the  nanopillar  due 
to  the  SPP-BW  mode,  while  ~  10%  is  absorbed  due  to  the  LSPR  mode  at  1295  nm.  It  should  be 
noted  that  the  absorption  efficiency  of  the  SPP-BW  mode  is  not  limited  to  18%  and  can  be 
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further  enhanced  by  optimizing  the  geometry  of  the  self  aligned  nanohole  lattice.  Fig.  1 1(b)  and 
Fig.  11(c)  shows  the  power  absorption  profile  of  the  SPP-BW  and  the  LSPR  modes  excited  in 
the  structure  at  normal  incidence.  In  the  case  of  the  SPP-BW  mode,  light  is  focused  away  from 
the  metal  into  the  nanopillar,  and  extends  downwards  from  the  gold  on  top  of  the  nanopillar.  In 
contrast  the  second  peak  at  1300  nm  has  an  absorption  profile  very  similar  to  the  LSPR  described 
in  Fig.  8(b),  which  demonstrates  that  LSPRs  supported  by  the  NOA  structure  can  be  utilized  for 
enhanced  absorption  in  an  array  format. 


Wavelength  [nm] 

Figure  14:  a)  Measured  angular  photoresponse  spectra  of  PENPD  showing  enhanced  photocurrent  due  to  both  SPP-1,0  and  LSPR  modes.  P 
=  580  nm  w  =  1 60  nm  h  =  320  nm  b)  Fitting  of  the  SPP-BW  and  LSPR  before  and  after  overlap  at  0“  and  30  0 


Two  peaks  can  be  observed  in  the  normal  incidence  photoresponse  spectrum.  The  first  peak  at 
923  nm  is  attributed  to  the  SPP-BW(0,i)  mode  corresponding  to  the  FDTD  mode  calculated  at 
1000  nm  while  the  second  peak  measured  at  1200  nm  is  attributed  to  the  LSPR  calculated  at 


1275  nm.  For  the  angular  photoresponse  characterization  the  sample  is  mounted  on  a  rotation 
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stage  described  previously  56  and  tilted  with  respect  to  the  incident  light.  When  the  device  is 
tilted  (0)  with  respect  to  nonnal,  incident  photons  have  an  in  plane  momentum  given  by 
kin-picme  =  k0  sin(6>)  and  the  phase  matching  condition  with  surface  plasmon  polaritons  is  altered. 
Fig.  5  shows  that  the  photoresponse  due  to  the  SPP-BW  mode  red  shifts  from  915  nm  at  0  =  0  0 
to  990  nm  at  0  =  30  0  and  spectrally  overlaps  with  LSPR  at  1180  nm  for  maximum  enhanced 
photoresponse.  The  peaks  of  the  photoresponse  spectra  are  fitted  to  Gaussian  peaks  to  detennine 
the  peak  position  and  line  width  of  the  SPP-BW  mode  and  LSPR  as  shown  in  Fig.  5(b).  A 
Gaussian  fitting  was  chosen  because  the  device  consists  of  ~  100,000  nanopillars  such  that  the 
photocurrent  measurement  is  an  ensemble  measurement.  At  zero  bias  0  =  0°  before  the  overlap 
of  the  SPP-BW  and  LSPR  modes  the  responsivity  at  1180  nm  due  to  the  LSPR  is  0.048  A/W. 
Due  to  the  overlap  of  the  SPP-BW  and  LSPRs  at  0  =  30  °  a  responsivity  of  0.096  A/W  is 
measured  at  1180  nm,  which  is  a  ~  x2  enhancement  in  the  responsivity  at  LSPR  wavelength  due 
hybridization  with  the  SPP-BW  mode. 

Plasmonic  bandstructure 

To  further  understand  the  coupling  of  the  hybrid  plasmonic  modes  excited,  photonic 
bandstructure  calculations  were  performed.  The  photonic  band  structure  of  the  device  is 
computed  by  the  fast  fourier  transfonn  (FFT)  of  the  time  decay  of  the  electric  field  intensity  time 
monitor  in  the  nanopillar.  An  excitation  dipole  placed  above  the  elongated  hole  with  Bloch 
boundary  conditions,  and  the  in-plane  momentum  is  varied  in  the  y-direction.  These  simulations 
performed  with  periodic  boundary  conditions  reveal  several  modes  which  did  not  previously 
appear  in  the  simulations  performed  with  PML  boundary  conditions  (Fig.  2(a)  )  -  we  consider 
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these  to  be  SPP-BW  modes  supported  by  the  gold/bcb  interface  and  the  gold/air  interface. 
However,  only  the  SPP-BW  mode  due  to  the  gold/BCB  interface  results  in  useful  optical 
absorption  in  the  nanopillar. 

Fig.  14  (a)  shows  the  FDTD  simulated  bandstructure  of  the  device  with  P  =  580  nm  ,  w  =  160 
nm,  and  h  =  320  nm.  At  the  gamma  point  light  couples  predominantly  into  a  band  with  a  very 
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Figure  15:  a)  Simulated  photonic  bandstructure  of  NP-PEPD  with  P  =  580  nm  showing  the  SPP-BW  Gold-BCB  modes  and  Gold-Air 
modes,  b)  the  absorbed  power  in  the  nanopillar. 

low  group  velocity  at  1000  nm  (~310  THz)  such  that  the  light  is  efficiently  absorbed  in  the 
nanopillar,  which  is  ideal  for  vertically  incident  photodetectors  60.  There  is  also  a  secondary  peak 
at  ~  1300  nm  (~  230  THz)  due  to  the  LSPR,  which  shows  little  dispersion  with  in  plane 
momentum  characteristic  of  a  localized  resonance.  Fig.  14(a)  also  shows  the  analytical  SPP-BW 
dispersion  of  a  gold  hole  array  on  polymer  overlaid  with  the  FDTD  simulated  bandstructure.  The 
analytical  expression  for  the  dispersion  of  the  SPP-BW  modes  is  given  by  Eq.  (5)  ’  ,  where 

^spp-bw,  's  the  resonant  scattering  orders  of  the  metal  hole  array,  P  is  the  pitch,  sd  is  the 
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dielectric  constant  of  the  polymer,  sAu  is  the  dielectric  constant  of  the  gold,  and  kinpkme  is  the  in¬ 
plane  momentum  in  the  direction  perpendicular  to  long  side  of  the  hole.  There  are  three  white 
lines  which  denote  the  Gold/polymer  SPP-BW  modes. 
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The  SPP(i,o)  and  SPP(.i,o)  modes  at  344THz  are  degenerate  at  the  gamma  point  and  does 
not  exhibit  in-plane  dispersion  as  these  modes  propagate  in  the  x-direction.  The  SPP(o,i)  mode  is 
also  at  344  THz  at  the  gamma  point  and  shifts  to  shorter  wavelengths  with  increasing  in  plane 
momentum.  However,  the  SPP(o,-i)  mode  shifts  towards  longer  wavelength  with  increasing  in 
plane  momentum,  and  is  the  candidate  SPP-BW  mode  for  overlapping  with  the  NOA  resonance 
for  enhanced  photoresponse.  Away  from  the  gamma  point  the  dispersion  of  the  analytically 
computed  SPP(o.-i)  matches  that  of  the  numerically  computed  bandstructure  as  is  seen  in  the 
dotted  white  line  in  Fig.  14(a).  However  this  prediction  fails  close  to  the  Gamma  point  (0  <  10 
degrees). 

The  dashed  red  lines  in  Fig.  14(a)  correspond  to  the  incident  illumination  being  tilted  an 
angle  9  with  respect  to  normal.  From  0°  to  20°  the  dispersion  of  the  FDTD  computed  electric 
field  intensity  at  301  THz  does  not  considerably  red-shift  with  in-plane  momentum,  and  is 
possibly  do  to  SPP-BW  mode  having  a  very  low  group  velocity.  However,  beyond  20°  the 
SPP-BW  red-shifts  further  such  that  at  30°  the  SPP-BW  has  the  greatest  overlap  with  LSPR.  The 
dispersion  of  the  numerically  computed  bandstructure  matches  the  analytical  expression  for  the 
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SPP-BWV-i)  mode.  Further  increasing  the  incident  angle  results  in  the  electric  field  intensity  due 
to  the  SPP(o,-i)  being  much  weaker  in  the  nanopillar  resulting  in  reduced  photoresponse.  The 
bandstructure  calculations  qualitatively  corresponds  to  the  angular  photoresponse  measurements 
shown  in  Fig.  5.  Due  to  the  low  group  velocity  of  the  SPP-BW  mode  the  photoresponse  peak 
does  not  shift  with  in  plane  momentum  upto  a  tilt  of  20°,  beyond  20°  the  SPP-BW  red  shifts.  In 
addition  to  the  Gold/BCB  SPP-BW  modes  the  structure  also  supports  SPP-BW  modes  due  to  the 
Gold/Air  interface  shown  as  the  dotted  gray  line,  which  has  a  good  fit  with  the  numerically 
computed  FDTD  bandstructure.  However,  these  modes  don't  contribute  to  electric  field  intensity 
enhancement  in  the  nanopillar  resulting  in  useful  photogeneration. 

Fig.  6(b)  shows  the  power  absorbed  in  the  nanopillar  resulting  in  useful  photoresponse. 
The  SPP-BW  mode  from  the  gold/air  interface  (gray  line  in  Fig.  6(a)  doesn’t  contribute  to  useful 
photogeneration.  However,  the  coupling  of  the  SPP-BW  (0,-1)  from  air/BCB  interface  and  the 
LSP  resonance  results  in  enhanced  photoresponse.  The  calculated  power  absorbed  spectra  in  Fig. 
6(b)  supports  the  experimental  result  that  a  combination  of  SPP-BW(o,-n  and  the  LSP  resonance 
is  producing  enhanced  photoresponse.  Hybridizing  SPP-BW  modes  with  LSPRs  results  in 
photoresponse  spectra  that  is  more  broadband  than  the  individual  plasmon  resonances. 
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3.0  GaAs  homojunction  3D-NOAADs 

Efficient  carrier  punch-through  from  the  low  field  absorption  region  into  the  high  field 
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multiplication  region  is  of  paramount  importance  in  the  design  of  the  3D  NOAADs. 
Photogenerated  carriers  can  recombine  due  to  surface  states  of  the  nanopillar  or  trap  states  at  the 
InGaAs-GaAs  heterojunction  limiting  device  efficiency.  Three-dimensional  finite  element 
modeling  of  the  nanopillar  heterojunctions  and  its  surfaces  are  a  necessary  component  in 


Figure  16:  Plasmonically  enhanced  photogeneration  rates  in  the  low  field  GaAs  absorber  a)  on  the  surface  b)  a  cut  through  the 
nanopillar 

nanostructure  device  design.  Three-dimensional  Drift-diffusion  to  model  3D  NOAAD  device 

characteristics,  by  solving  the  basic  Poisson  and  drift-diffusion  equations  without  simplification 

in  three  dimensions  and  can  be  used  for  arbitrary  device  geometries.  New  physical  models  based 

on  specific  material  properties  of  the  nanopillar  can  be  incorporated  such  as  trap  assisted 

recombination  and  tunneling.  Figure  15  (a)  and  (b)  shows  electron  carrier  density  due  to  the 

plasmonically  enhanced  photogeneration.  Carrier  extraction  might  be  reduced  due  to  electron 

trapping  at  the  InGaAs  -  GaAs  heterojunction,  we  believe  this  can  be  overcome  with  a 
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compositionally  graded  InGaAs  absorption  region  reducing  the  band  offset  for  carrier  extraction. 
The  3D  photogeneration  rates  computed  by  electromagnetic  simulation  will  be  imported  into  the 
3D  drift  diffusion  calculations  to  estimate  bias  dependent  quantum  efficiency  of  the  3D 
NOAAD.  Models  of  dark  current,  trap  densities  and  estimate  carrier  recombination  at  the 
surface  will  be  extracted  by  fitting  to  temperature  dependent  dark  current  measurements  to 
experimental  results.  The  3D  electric  field  profiles  generated  by  drift  diffusion  modelling  will 
also  used  to  model  avalanche  gain  and  excess  noise  of  the  3D  NOAADs. 

3.1  GaAs  core -shell  nanopillar  APDs 

Figure  16a)  illustrates  the  structure  of  the  NP  APD  comprised  of  a  p-GaAs  absorber  (150 
nm  in  width  and  900  nrn  in  height)  atop  a  core-shell  GaAs  PN  junction.  The  n-GaAs  core  is  80 
mn  in  width  and  800  nm  in  height,  while  the  p-GaAs  shell  is  35  nm  thick  surrounding  the  core. 

1  Q  _ O 

The  intended  n-core  and  p-shell  doping  levels  are  both  ~  10  cm  . 
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Figure  17:  a)  Schematic  of  the  separate  absorption  multiplication  avalanche  diode,  (b)  Self-aligned  metal  nanohole  lattice  for 
plasmonically  enhanced  absorption,  (c)  Plasmonicallv  enhanced  photogeneration  within  the  nanopillar. 


The  nanopillars  are  encapsulated  with  an  kio.5Gao.5P  high  bandgap  passivation  layer  (not 

fx'X 

shown)  to  mitigate  the  effect  of  surface  states  .  Figure  16  b)  shows  a  SEM  image  of  the  3D  self 
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aligned  plasmonic  antenna  for  enhanced  absorption  within  the  GaAs  nanopillar,  engineered  to 
take  place  within  the  GaAs  bandgap  at  X  ~  730  mn.  Photogenerated  electrons  created  within  the 
tip  of  the  nanopillar  drift  to  the  core-shell  multiplication  region  with  applied  bias  as  shown  in 
Fig.  16  c). 

Figure  17a  shows  the  typical  light  and  dark  current-voltage  ( I—V)  characteristics  of  the 
nanopillar  array  contacting  -149160  NPs  at  300  K.  The  dark  I-V  characteristic  in  Figure  17a  can 
be  fit  to  the  ideal  diode  equation  where  the  current,  /,  is  given  by  1  =  ID  exp(g(V  -  IRs)/nkT),  ID  is 
the  saturation  current,  Rs  is  the  series  resistance,  q  is  the  elementary  charge,  T  is  the 
temperature,  V  is  the  bias  voltage,  and  ideality  factor  n.  The  nanopillar  array  devices  exhibit 

excellent  diode  characteristics  with  a  low  series  resistance  of  -40  f2,  ideality  factor  of  1.6,  and 

—12 

saturation  current  of  3  x  10  A.  These  diode  characteristics  suggest  a  high  quality  PN  junction 
is  fonned  within  the  nanopillar.  At  zero  bias,  under  light  conditions,  an  extremely  low 
photocurrent  of  20  nA  is  measured,  indicating  that  a  short  diffusion  length  ( «  1  pm)  results  in 
the  photogenerated  carriers  not  being  extracted  by  the  core-shell  electric  field.  However,  as  the 
reverse  bias  is  increased,  the  photocurrent  increases  corresponding  to  both  (1)  an  increase  in  the 
field  dependent  diffusion  length  of  photogenerated  carriers64  and  (2)  an  extension  of  the 
depletion  region  into  the  low  field  absorber.  At  -3  V,  all  of  the  photogenerated  carriers  in  the 
nanopillar  are  extracted,  forming  a  plateau  region  in  the  photocurrent  resulting  in  a  responsivity 
of  0.27  A/W  at  730  mn.  This  bias  voltage  corresponds  to  the  punch  through  voltage  indicating 
the  maximum  expansion  of  the  depletion  region.  Further  biasing  the  device  results  in  the 
photocurrent  being  multiplied  beyond  the  primary  photocurrent  due  to  avalanche  multiplication. 
Figure  2b  shows  the  gain  characteristic  calculated  by  subtracting  the  dark  current  from  the  light 
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current  and  normalizing  with  respect  to  the  primary  photocurrent  at  the  onset  of  the  plateau 
region.  The  unity  gain  primary  photocurrent  is  measured  at  -3.2  V.  Increased  reverse  bias  results 
in  a  mono  tonic  increase  in  avalanche  gain  until  gains  as  high  as  -216  can  be  achieved  at  -12  V 
due  to  avalanche  multiplication  in  Fig.  17b). 


Figure  18:  a)  Light  and  dark  I-V  characteristic  and  b)  the  multiplication  of  GaAs  core-shell  diode 

A  deeper  understanding  of  the  measured  gain  characteristics  and  a  correlation  with  the 
volume  of  the  avalanche  region  can  be  obtained  from  modeling  the  capacitance-voltage  (C— V) 
measurements  using  the  doping  and  the  geometry  of  the  nanopillar.  The  best  fit  for  the  C—V 
characteristic  can  be  used  to  estimate  the  high  electric  field  volume  responsible  for  the  avalanche 
gain.  Figure  18  shows  capacitance  (C—V)  characteristics  as  a  function  of  area  at  a  high  frequency 
of  1  MHz  to  sample  purely  the  junction  capacitance.  The  capacitance  scales  with  area  indicating 
that  the  junction  capacitance  of  each  nanopillar  is  acting  in  parallel  and  a  very  high  uniformity  in 
the  formation  of  the  core-shell  PN  junction  in  the  nanopillar  array.  At  zero  bias  the  capacitance 
due  to  the  core-shell  junction  is  -0.347  fF/prn  .  A  small  reverse  bias  of-1  V  results  in  a  -50% 
drop  in  the  capacitance,  until  higher  reverse  biases  results  in  the  saturation  of  the  capacitance 
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to  -0.05  fF/pm  '  We  roughly  estimate  that  this  ultrasmall  saturation  capacitance  would  allow 
nanopillar  diode  arrays,  with  as  few  as  3000  NPs  to  have  a  RC  limited  bandwidth  exceeding  40 
GHz. 


Figure  19:  a)  Capacitance-Voltage  (CV)  characteristics  of  nanopillar  arrays  with  different  number  of  nanopillars  b)  Extracted 
average  capacitance  per  nanopillar  compared  with  a  modelled  capacitance  c)  Schematic  of  the  structure  used  for  3D  modeling  of 
capacitance 

Figure  18(b)  shows  the  average  capacitance  per  nanopillar  overlaid  with  C-V  modeling 
(Curves  1  and  2)  of  the  NP  core-shell  junction.  The  average  junction  capacitance  per  nanopillar 
can  be  estimated  by  subtracting  the  contribution  to  the  capacitance  due  to  BCB  and  dividing  by 
the  number  of  NPs  in  the  array.  The  modeled  C-V  curves  are  generated  by  solving  the  3D 
Poisson  equations  65  taking  into  account  the  core-shell  geometry  and  doping  of  the  hexagonal 
nanopillar.  In  both  Curves  1  and  2  the  modeled  core  is  80  nm  in  diameter  and  800  nm  in  height 
and  uniformly  doped  n-type  at  a  doping  concentration  of  3x10  cm'  .  The  doping  of  the  p-type 
doped  shell  is  coarsely  graded,  i.e.  divided  into  three  sections  including  base,  middle  and  tip 
accounting  for  expected  growth  rate  changes  during  nanopillar  epitaxy.  ’  The  doping 

173173 

concentration  within  the  base,  middle  and  tip  sections  are  constant  at  3x10  cm'  ,7x10  cm' 
and  1x10  cm'  ,  respectively.  However,  the  CV  characteristics  are  very  sensitive  to  the 
thickness  of  the  middle  section  surrounding  the  tip  of  the  core.  Thus,  Curve  1  (2)  represent  a 
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base  section  spanning  500  mn  (600  nm),  the  middle  section  of  150  nm  (100  nm),  and  the  tip  at 
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Figure  20:  a)  Contour  plot  of  the  electric  field  in  the  nanopillar  at  -10  V  b)  Axial  electric  field  profile  along  the  pillar 
as  a  function  of  reverse  bias,  c)  Calculated  gain  using  local  and  DSMT  compared  with  experimental  gain  d)  Shows  the 
characteristic  temperature  dependence  of  avalanche  breakdown 

150  nm  (100  nm). 


A  clear  trend  can  be  seen  in  that  the  peak  capacitance  at  low  bias  voltages  drops  as  the  size  of 
the  high  doped  region  encapsulating  the  tip  of  the  core  is  reduced,  while  the  saturation 
capacitance  also  increases.  Curve  1  is  chosen  to  most  accurately  represent  the  characteristics  of 
the  device  in  terms  of  accurately  predicting  the  punch  through  voltage,  the  shape  of  the  CV 
characteristic  and  the  saturation  capacitance  Csat  ~  9  aF.  The  CV  characteristic  of  Curve  1  can  be 
understood  as  the  radial  expansion  of  the  3D  depletion  region  by  initially  depleting  the  core-shell 
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junction  and  then  extending  axially  to  the  pillar  edge.  At  the  bias  voltage  of  -5  V  the  depletion 
region  has  expanded  axially  resulting  in  the  saturation  of  the  capacitance  at  -  9  aF  for  a  single 
nanopillar. 

Figure  19(a)  shows  a  contour  plot  of  the  electric  field  through  the  center  of  the  nanopillar 
(in  the  x-plane)  at  the  breakdown  voltage  of  -10V.  The  graded  doping  of  the  p-shell  results  in 
the  high  field  (E  >  5xl05  V/cm)  region  being  confined  to  within  310  run  of  the  tip  of  the  core. 
Since  very  high  electric  fields  are  confined  to  within  a  very  small  region  of  the  core-shell 
junction  the  resulting  impact  ionization  is  also  being  confined  to  a  very  small  volume.  The  thin 
avalanche  volume  results  in  the  measured  gain  characteristics  to  be  dominated  by  dead  space 
effects. 

Figure  19(b)  shows  the  axial  electric-field  profiles  as  a  function  of  reverse  bias  through 
the  center  of  the  nanopillar  (at  the  location  of  the  dotted  white  line  in  Figure  19(a))  for  the 
structure  described  by  Curve  1 .  The  avalanche  region  at  the  breakdown  voltage  is  confined  to  an 
axial  region  -310  nrn  in  width  with  highly  localized  fields  as  high  as  1  MV/cm  at  the  tip  of  the 
nanopillar  core.  The  bias  dependent  electric  field  profiles  can  be  used  to  model  the  expected  gain 

rn  ...  Q  . 

using  both  local  theory  for  carrier  multiplication  and  dead  space  multiplication  theory  .  It  is 
expected  that  modeling  the  gain  in  terms  of  dead  space  multiplication  theory  (DSMT)  would 
lower  the  mean  gain  at  a  given  bias  and  shift  the  breakdown  voltage  to  higher  biases.  In  both 
cases,  it  is  assumed  that  a  majority  of  carriers  contributing  to  avalanche  multiplication  traverse  a 
path  through  the  center  of  the  nanopillar  as  shown  by  the  white  dotted  line  in  Figure  19(a).  In 
terms  of  the  local  theory  the  mean  multiplication  ( M  )  is  given  using  Eq.  1  below,  where  OC  is 
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the  electron  ionization  coefficient,  [3  is  the  hole  ionization  coefficient,  and  is  evaluated  at 
x  =  W  where  W  is  the  position  of  the  depletion  edge. 
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For  the  case  of  DSMT  the  electric-field  profiles  in  Fig.  19(b)  were  also  used  to  calculate 
the  dead-space  profiles  across  the  nanopillar  and  the  non-local  ionization  coefficients  for 
electrons  and  holes  in  GaAs  as  described  by  Saleh  et  al. 14  With  this  information  we  calculate  the 
probability  density  functions  of  the  random  distance  of  ionization  for  electrons  and  holes.  A 
hard-threshold  dead  space  is  assumed,  for  which  the  ionization  coefficient  is  zero  for  carriers  that 
have  traveled  a  distance  shorter  than  the  dead  space  and  is  a  function  of  the  electric  field  for 
carriers  that  have  traveled  a  longer  distance. 

Figure  19(c)  shows  the  calculated  gain  using  local  theory  and  DSMT  using  the  electric 
fields  of  Curve  1  overlaid  with  the  experimental  gain.  Local  theory  is  shown  to  produce  an  over 
estimation  to  the  measured  gain  resulting  in  an  under  estimation  of  the  breakdown  voltage. 
However,  the  mean  gain  estimated  by  DSMT  is  able  to  predict  the  experimental  gain  and 
breakdown  voltage  more  accurately  than  the  local  theory,  since  the  effect  of  the  dead  space  are 
significant. 

Figure  19(d)  shows  the  temperature  dependent  I-V  characteristics  of  the  core-shell  NP 

diodes  at  high  reverse  biases  from  77K  to  300K.  All  curves  show  similar  behavior:  a  clear 

transition  from  a  generation-recombination  dominated  current  to  an  avalanche  current  close  to 

the  breakdown  voltage  at  ~  9.5  V.  The  breakdown  voltage  is  consistent  with  the  predictions  of 
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the  DSMT.  This  dramatic  change  in  the  I-V  characteristic  suggests  that  dark  carriers  injected  into 
the  high-field  core-shell  region  are  multiplied.  The  breakdown  voltage  can  be  calculated  by 
extracting  the  voltage  at  which  multiplication  (M)  goes  to  infinity.  The  unmultiplied  dark  current 
is  taken  to  be  -7  V  from  which  the  multiplication  factor  M  is  calculated.  The  intersection  of  1/M 
vs.  bias  is  used  to  estimate  the  breakdown  voltage  at  each  temperature  (See  Supporting 
Infonnation).  The  breakdown  voltage  shifts  to  larger  biases  with  increasing  temperature  with  a 
shift  of  +3.3  mV/K,  such  a  small  shift  in  the  breakdown  voltage  suggests  that  the  avalanche  gain 
is  produced  by  a  high  field  region  on  the  order  of  the  ionization  path  length  of  electrons  as  higher 
temperatures  result  in  greater  energy  scattering  of  carriers.  The  clear  shift  in  the  breakdown 

voltage  is  in  stark  contrast  to  previous  reports  of  GaAs  nanowire  APDs  where  a  positive  shift  in 

68 

the  breakdown  voltage  is  not  clearly  observed. 

Excess  noise  characterization  of  GaAs  3D-NOAADs 

GaAs  3D  NOAADs  were  characterized  for  gain  and  excess  noise  factors  in  collaboration 
with  Prof.  Siddhartha  Ghosh  at  University  of  Illinois  Chicago.  Devices  were  probed  using  a  Janis 
500  T  cryogenic  probe  station  at  300  K  such  that  DC  bias  was  applied  using  a  bias-tee  and  AC 
noise  was  measured  using  a  Agilent  8973b  noise  meter.  Optical  illumination  was  provided  using 
a  He-Ne  laser  calibrated  to  output  ~  200  pW.  Unity  gain  was  calculated  from  Fig.  13b)  after 
punch  through  of  the  depletion  region  into  the  absorption  region. 
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Figure  20  shows  the  measured  excess  noise  factor  a  function  of  multiplication  gain  in  the 
GaAs  nanopillar  avalanche  detectors.  At  low  gains  the  extracted  excess  noise  factor  is  dominated  by 


Figure  21:  Excess  noise  factor  vs.  Gain  for  GaAs  NOAAD  fitted  with  Mcyntyre  fitting  of  keff 

the  output  noise  of  the  preamplifier.  However,  at  higher  gains  the  excess  noise  factor  approaches  fits 
to  the  McIntyre  equation.  Several  devices  from  the  same  nanopillar  growth  run  have  been  measured 
and  similar  low  excess  noise  has  been  observed.  The  measured  excess  noise  factors  are  well  below 
the  ratio  of  the  ionization  coefficients  reported  for  GaAs  bulk  avalanche  detectors.  The  effective  ratio 
of  ionization  coefficients  for  electrons  to  holes  is  measured  to  be  less  than  0.1.  For  bulk  GaAs,  the 
ratio  is  close  to  1\  Low  excess  noise  factors  (keff  ~  0.3)  have  been  measured  in  bulk  GaAs  diodes 
where  the  multiplication  region  ~  100  nm.  The  physical  origin  of  such  a  low  excess  noise  factor  is 
still  being  investigated  using  dead  space  multiplication  theory  (DSMT). 
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4.0  InGaAs-GaAs  heterojunction  NOAADs 

Incorporating  indium  (In)  within  the  nanopillar  during  the  GaAs  nanopillar  allows  the 
extension  of  the  cut-off  wavelength  to  the  1-1.3  pm  range41.  Lattice  mismatched  InGaAs 
nanopillar  arrays  were  grown  by  selective  area  epitaxy  metal  organic  vapor  phase  epitaxy  on 
GaAs  substrates  as  shown  in  Fig.  21  a).  The  nanopillars  were  passivated  by  wide  bandgap 
InGaP  shells  to  reduce  surface  related  dark  current  densities  resulting  in  unifonn  arrays  as  shown 
in  Fig.  21  b)  and  c).  The  nanopillar  arrays  were  fabricated  into  high  speed  APDs  with  GSG 
contact  pads  utilizing  self  aligned  metal  hole  array  antenna  for  the  signal  contact  as  shown  in 
Fig.  2 Id)  and  e). 


Figure  22:  a)  Schematic  of  nanopillar  structure  investigated  for  1  micron  nanopillar  APD  b)  SEM  of  10  pm  x  10  pm  array  c)  Close 
up  view  of  InGaAs  nanopillar  passivated  with  InGaP  d)  Close  up  SEM  of  final  device  showing  plasmonic  metal  hole  array  antenna 
e)  Nanopillar  array  integrated  with  GSG  contact  pads 

properties  of  the  APDs  were  investigated  by  current  voltage  (IV)  measurements  under  dark  and 

under  light  as  shown  in  Figure  22.  Sample  AF024  was  grown  such  that  the  Indium  composition 
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1 7  7 

36%,  p-type  Zinc  doping  of  ~  10  cm'  .  Sample  AF025  has  an  Indium  composition  36%,  p-type 
Zinc  doping  of  ~  10  cm'  .  Sample  AF026  has  a  Indium  composition  27%,  p-type  Zinc  doping  of 
~  1017cm'3 

All  photocurrent  measurements  perfonned  with  a  fiber  coupled  1060nm  laser  calibrated 
to  21  |JW  of  incident  optical  power.  The  device  characteristics  of  three  growths  are  reported 
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Figure  23:  Light  dark  current  voltage  and  gain  characteristics  of  InGaAs  nanopillar  APDs 


here.  The  dark  IV  characteristics  show  that  increasing  the  doping  of  the  nanopillar  lowers  the 
breakdown  voltage  (~  9  V  )  with  increased  dark  current  (2.38  nA/pillar  @  M  =  10),  while  lower 
doping  concentrations  produces  higher  breakdown  voltages  (~  20  V)  with  lower  dark  currents 
(1.37  pA/pillar  @  M  =  10).  Reducing  the  indium  composition  of  the  nanopillar  is  also  shown  to 
reduce  dark  current  (143  fA/pillar  @  M  =  10)  and  increase  the  breakdown  voltage.  In  terms  of 
the  performance  of  the  devices  under  light,  the  best  combination  of  gain  and  quantum  efficiency 
was  observed  for  device  in  Fig.  2b)  with  a  unity  gain  quantum  efficiency  of  21%  and  a 
maximum  gain  of  1430. 

4.1  Electric  field  of  InGaAs -NO AADs 

As  shown  earlier  the  avalanche  breakdown  characteristics  of  InGaAs  NOAADs  are 

heavily  dependent  of  the  doping  profile  of  the  nanopillar.  A  detailed  analysis  of  the  capacitance 
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voltage  characteristics  are  performed  to  extract  doping  concentrations  within  the  nanopillar.  The 
InGaAs  nanopillar,  GaAs  substrate,  the  geometry  of  the  top  contact  and  the  in-filling  BCB  layer 
is  modeled  in  Sentaurus.  The  doping  concentration  of  the  InGaAs  nanopillar  is  varied  from  1  x 
10  cm'  at  the  base  to  5  x  10  cm'  at  the  tip.  The  extracted  fit  to  the  experimental  capacitance 
is  shown  in  Fig.  20  a).  Fig.  20b)  shows  the  3D  electric  field  distribution  of  the  InGaAs  nanopillar 
GaAs  substrate  heterojunction.  It  should  be  noted  that  the  BCB  planarization  layer  smears  out 
the  field  at  the  substrate  nanopillar  junction,  similar  to  the  breakdown  characteristics  observed  by 
in  silicon  nanowire  pn  junctions1.  Since  the  p-doped  InGaAs  region  has  a  doping  gradient  a  quasi 
electric  field  is  created  which  aids  in  extraction  of  photogenerated  electrons  into  the 
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Figure  25:  a)  Fitting  of  the  experimental  CV  characteristics  b)  3D  electric  field  profile  including  dielectric  polymer  layer  c)  1-D 
electric  through  the  center  of  the  nanopillar.  Inset  shows  a  electric  field  in  the  x-z  plane 

multiplication  region.  At  8  V  reverse  bias  the  high  field  region  is  confined  to  within  a  ~  200  mn 
from  the  base  of  the  pillar  as  shown  in  Fig.  22c).  The  drift  electric  field  in  the  absorption  region 
reaches  ~  400  V/cm,  which  is  sufficiently  high  to  cause  carriers  to  drift64  into  the  multiplication 
region,  instead  of  recombining. 
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4. 1  Excess  noise  characterization  oflnGaAs  NOAADs 

InGaAs  NOAADs  are  characterized  for  excess  noise  gain  characteristics  using  custom  setup  built 
at  UCLA.  The  schematic  of  the  excess  noise  setup  is  shown  in  Fig  XX  below.  The  main 
technical  challenge  to  be  overcome  in  making  noise  measurements  is  to  accurately  attribute  the 
measured  noise  to  the  device  under  test  as  opposed  to  external  noise  sources  such  as  DC  power 
supplies,  preamplifier  noise,  unshielded  cables  and  laser  induce  Relative  Intensity  Noise.  High 
speed  GSG  contact  pads  reduces  the  contribution  to  the  noise  signal  from  external  sources.  Low 
noise  DC  bias  is  provided  by  the  Agilent  B2962A  through  the  picosecond  pulse  labs  model  5543 
bias  T.  The  AC  current  noise  is  amplified  by  the  low  noise  current  preamplifier  Stanford 
Research  model  SR570  and  spectral  noise  density  is  displayed  using  the  Stanford  Research  FFT 
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Figure  26:  Noise  Spectral  Density  (NSD)  vs.  photocurrent  b)  Excess  noise  factor  vs.  Gain  (M) 

signal  analyzer  model  SR780.  The  noise  current  density  is  extracted  from  the  "flat  region"  of  the 
noise  spectrum  beyond  the  1/f  noise  knee.  Light  is  coupled  into  the  APD  using  a  multimode 
lensed  fiber  with  a  spot  size  ~  50  pm  coupled  to  a  pigtailed  fiber  laser  with  a  FC  termination. 
Shot  noise  of  the  DC  photocurrent  is  confirmed  below  unity  gain  as  shown  in  Fig.  23  a)  fitting  the 
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relationship  im  =  ^2 ql ph  .  With  increasing  photocurrent  the  current  noise  spectral  density 

increases  beyond  the  shot  noise  floor.  The  increased  noised  can  be  quantified  in  tenns  of  an 

5 

excess  noise  factor  (F),  such  that  F  = - — — - 69 .  The  onset  of  avalanche  gain  at  low  electric 

2  qIphM2 

fields  can  be  hidden  by  the  movement  of  the  depletion  edge  with  increased  bias  improving 
carrier  collection.  The  bias  voltage  for  the  additional  noise  beyond  the  shot  noise  floor  was  used 
to  identify  unity  gain  .  Excess  noise  factors  for  the  InGaAs  3D  NOAADs  are  significantly  less 
that  bulk  APDs  where  keff  approached  unity  as  the  multiplication  region  width  is  less  that  1 
micron. 

4.  2  Gain-bandwidth  characterization  of  InGaAs  NOAADs 

The  bandwidth  of  APDs  can  either  be  transit-time  limited  or  RC-limited  depending  on 

the  fabrication  scheme  of  the  device.  Emmons71  showed  for  transit  time  limited  APDs  the 

bandwidth  would  reduce  as  soon  as  the  gain  (M)  >  1  /  keff.  However,  in  the  case  of  the  InGaAs 

NOAADs  the  devices  were  not  fabricated  on  semi -insulating  substrates  to  reduce  parasitic 

capacitances,  and  therefore  RC-limited.  The  signal  contact  being  deposited  on  top  of  the  BCB 

creates  a  parasitic  capacitance  of  ~  280  fF,  together  with  a  series  resistance  of  ~  220  Q  produces 

a  RC  limited  bandwith  of  2.6  GHz.  Impulse  response  and  associated  DC  gain  of  typical  InGaAs 

NOAAD  devices  are  shown  in  Fig.  26  below.  DC  gains  as  high  as  982  are  achieved  at  a  reverse 

bias  of  7.6  V.  From  the  impulse  response  measurements  shown  in  Fig.  26  b)  the  response  of  the 

InGaAs  NOAADs  are  non-gaussian  with  a  fast  rise  time  (~17  ps)  and  a  significantly  slower  fall 

time  (~  90  ps).  The  bandwidth  of  the  InGaAs  NOAAD  is  essentially  limited  by  the  fall  time, 

which  is  attributed  to  the  RC  time  constant  of  parasitic  BCB  capacitance.  A  maximum  gain 
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bandwidth  product  of  2.8  GHz  x  711  =  2  THz.  Such  high  gain-bandwidth  products  is 

70 

exceptionally  high  for  the  InGaAs  material  system  . 


Figure  27:  Dark/light  IV  characteristics  of  InGaAs  NOAAD  b)  Impulse  response  of  InGaAs  NOAAD  to  80  fs  Ti-sapphire  laser 
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5.0  Future  outlook 

The  devices  realized  in  this  thesis  demonstrate  the  proof  of  concept  of  NOAADs  in 
accomplishing  wavelength  tunable,  low  operating  bias,  high  gain,  low  excess  noise  and  high 
speed  avalanche  photodetectors.  Further  improvements  can  be  made  to  engineer  the  electric  field 
of  the  NOAAD,  such  that  high  gains  can  be  achieved  with  lower  dark  currents.  The  mechanisms 
for  dark  current  generation  needs  to  be  studied  by  varying  the  nanohole  diameter  to  isolate  bulk 
and  surface  induced  dark  current.  This  is  particularly  challenging  in  nanopillars  grown  by 
selective  area  epitaxy  since  the  growth  conditions  of  the  nanopillar  changes  with  the  size  of  the 
electron  beam  defined  nanohole.  InP  and  InAlAs  nanopillars  can  also  be  explored  as  candidate 
materials  for  multiplication  regions,  which  would  potentially  reduce  dark  currents.  In  terms  of 
commercial  application  to  LADAR  focal  plane  arrays,  NOAADs  need  to  be  integrated  with 
planar  epitaxy  on  (111)  substrates  to  achieve  hybridization  with  silicon  Read  out  Integrated 
Circuits. 
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